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Abstract. This paper explores the evolution of three-dimensional ob- 
jects with a simple generative encoding, known as the Superformula. 
Evolving three-dimensional objects has long been of interest in a wide ar- 
ray of disciplines, from engineering (e.g., robotics) to biology (e.g., study- 
ing morphological evolution). While many representations have been 
presented, ranging from direct encodings to complex graphs and gram- 
mars, the vast majority have possessed complex underlying encodings, 
which were necessary to produce varied morphologies. Here we explore 
the target-based evolution of Superformula, which is simply encoded as 
a vector of reals and show that it is possible to generate very closely 
matching designs of a number of complex three-dimensional objects. 



1 Introduction 

The evolution of geometric models to design arbitrary three-dimensional mor- 
phologies has been widely explored. Early examples include Watabe and Okino's 
lattice deformation approach [28] and McGuire's sequences of polygonal opera- 
tors [T7] . Sims [IS] evolved the morphology and behaviour of virtual creatures 
that competed in simulated three-dimensional worlds with a directed graph en- 
coding. Bentley [3] investigated the creation of three-dimensional solid objects 
via the evolution of both fixed and variable length direct encodings. The ob- 
jects evolved included tables, heatsinks, penta-prisims, boat hulls, aerodynamic 
cars, as well as hospital department layouts. Eggenberger [6] evolved three- 
dimensional multicellular organisms with differential gene expression. Jacob and 
Nazir [15] evolved polyhedral objects with a set of functions to manipulate the de- 
signs by adding stellating effects, shrinking, truncating, and indenting polygonal 
shapes. More recently, Jacob and Hushlak [M] used an interactive evolutionary 
approach with L-systems [21j to create virtual sculptures and furniture designs. 

Few evolved designs however have been manufactured into physical objects. 
Conventionally evolved designs tend to be purely descriptive, specifying what to 
build but not how it should be built. Thus there is always an inherent risk of 
evolving interesting yet unbuildable objects. Moreover, high fidelity simulations 
are required to ensure that little difference is observed once the virtual design 
is physically manifested. In highly complex design domains such as dynamic 



objects, the difference between simulation and reality is too large to manufacture 
designs evolved under a simulator. 

Funes and Pollack !8^ performed one of the earliest attempts to physically 
instantiate evolved three-dimensional designs by placing physical LEGO bricks 
according to the schematics of the evolved individuals. A direct encoding of 
the physical locations of the bricks was used and the fitness was scored using 
a simulator which predicted the stability of the composed structures. Pollack 
et al. [IH] subsequently evolved robot morphologies under simulation using a 
direct encoding and L-systems. The objects were ultimately fabricated by hand, 
and later, a three-dimensional printer. Additionally, Hornby and Pollack [13) 
used L-systems to evolve furniture designs, which were then manufactured by 
a three-dimensional printer. They found the generative encoding of L-systems 
produced designs faster and with higher fitness than a non-generative system. 
Generative systems are known to produce more compact encodings of solutions 
and thereby greater scalability than direct approaches (e.g., see [2^). 

Lohn et al. [16] evolved and manufactured X-band satellite antenna for 
NASA's ST5 spacecraft, representing the world's first artificially evolved hard- 
ware in space. Significantly, the antenna's performance was similar to a hand- 
designed antenna produced by an antenna-contractor. A generative encoding L- 
system was used which drew the three-dimensional model of the antenna when 
interpreted. The fitness was determined through the use of an antenna simulator. 

Recently, the generative encoding Compositional Pattern Producing Net- 
works (CPPNs) [26] have been used to evolve three-dimensional objects which 
were ultimately fabricated on a three-dimensional printer jll2| [S] . Both interac- 
tive and target-based approaches were explored. 

Significantly, most approaches have used simulations to provide the fitness 
scores of the evolved designs. Embodied evolutionary computing has typically 
referred to the existence of a physical solution in the fitness evaluation, and 
can be traced back to the origins of the discipline: the first Evolution Strate- 
gies were used to design jet nozzles as a string of real-valued diameters, which 
were then machined and tested for fitness [35]. Recently, Rieffel and Sayles [25] 
evolved circular two-dimensional shapes where each design was fabricated on 
a three-dimensional printer before assigning fitness values. Interactive evolution 
was undertaken wherein the fitness for each printed shape was scored subjec- 
tively. Each individual's genotype consisted of twenty linear instructions which 
directed the printer to perform discrete movements and extrude the material. 
As a consequence of performing the fitness evaluation in the environment, i.e., 
after manufacture, the system as a whole can exhibit epigenetic traits, where 
phenotypic characteristics arise from the mechanics of assembly. One such ex- 
ample was found when selecting shapes that most closely resembled the letter 
'A'. In certain individuals, the cross of the pattern was produced from the print 
head dragging a thread of material as it moved between different print regions 
and was not explicitly instructed to do so by the genotype. 

In this paper, we explore the evolution of a simple generative encoding to pro- 
duce three-dimensional designs for manufacture on a three-dimensional printer. 



2 Superformula 



Gielis f9'10| found that the forms of a large variety of plants and other living 
organisms can be modelled by a single, simple, geometric equation, forming a 
generalisation of a hyper-ellipse, termed the Superformula. Modifying the set of 
real- valued parameters to the Superformula generates myriad and diverse natu- 
ral polygons with corresponding degrees of freedom. The Superformula can be 
used to create three-dimensional objects, Supershapes, using the spherical prod- 
uct of two superformulas; in fact, by multiplying additional superformulas it can 
be extended to A^-dimensions. "In general, one could think of the basic Super- 
formula as a transformation to fold or unfold a system of orthogonal coordinate 
axes like a fan. This creates a basic symmetry and metrics in which distances can 
further be deformed by local or global transformations. Such additional transfor- 
mations increase the plasticity of basic Supershapes" pU]. Gielis' Superformula 
can be further generalised to increase the degrees of freedom, adding twist and 
further rotations, permitting the creation of more complex three-dimensional 
forms, including shells, mobius strips, and umbilic tori. Gielis' Superformula, 
which defines a supershape in 2 dimensions is given in the following equation, 
where r is the radius; cj) is the angle; a > 0, 6 > control the size of the su- 
pershape and typically = 1; and m (symmetry number), ni, n2 and (shape 
coefficients) are the real- valued parameters: 

r - fW I ^ (1) 

Using the spherical product, the extension to three dimensions: 

X — ri{9) X cosiO) X ^2(1^9) x cos{ip) 

y = ri{9) X sin{6) x r2{(p) x cos{ip) 
z — r2{(p) X sin(ip) 

Where — 7r/2 < (p < 7t/2 for latitude and — tt < < tt for longitude. 

Example shapes generated with the Superformula can be seen in Figure [T] 
where the cube, star, and heart can be generated from the same set of eight 
real- valued parameters; the torus requiring two additional parameters; the shell 
a total of twelve; and the mobius strip a total of fifteen. 

A "Superdupershape explorer" and its source code, licensed under Creative 
Commons Attribution- Share Alike 3.0 and GNU GPL license, can be found at 
[http : / / openprocessing . org/visuals/?visualID=2638 

Given a target shape it is often very useful to identify a representative for- 
mula. Optimisation methods, such as the Levenberg-Marquardt (LM) theory 
PU] , have typically been used to identify the best fitting supershape parameters 
(e.g., [H]). However LM cannot retrieve all of the parameters required for su- 
pershape fitting. Bokhabrine et al. [4] used a Genetic Algorithm (GA) [12] to 
evolve all supershape parameters for surface reconstruction (i.e., a target-based 



approach) using an inside-outside function [7 for fitness computation. Voisin 
et al. [37] later extended this to utihse a pseudo-Euchdean distance for fitness 
determination, yielding improved performance. Additionally, Morales et al. |18) 
used a GA to evolve TV-dimensional Superformula for clustering. 

3 Target-Based Evolution of Supershapes 

The cube, star, and heart shapes (as seen in Figure [T]) are here converted into 
125,000 voxel binary three-dimensional grids and used as the desired targets, 
where the fitness of an individual is the fraction of voxels that correctly match. 
The genotype of each individual in the population consists of eight real-valued 
parameters in the range [0,20] which affect the Superformula, giving rise to the 
Supershape. The GA proceeds with a population of 200 individuals, a per allele 
mutation rate of 25%, and mutation step size of ±rand{5) , where rand selects a 
real- valued number in the range [0,5]; the GA tournament size for both selection 
and replacement is set to 3. 

Figure [2] shows the fraction of total voxels matched to the target shapes 
during evolution of the Supershapes; results presented are an average of 10 ex- 
periments. Similar to a large number of voxels are quickly matched, however 
here the target object is not identifiable until approximately 99% are set cor- 
rectly. As such, the small differences in fitness between the treatments represent 
substantial differences in whether the target object is recognisable. In all cases, 
greater than 99.5% fitness is achieved. From Figure it can be seen that, on 
average, the GA takes approximately 1100 evaluations to reach >99% matching 
voxels of a target cube object and 3700 evaluations to achieve >99.9%. Fig- 
ure [2b] shows that on average approximately 3900 evaluations are required to 
reach >99% matching voxels of a target star object and 16100 evaluations to 
achieve >99.5%. Finally, Figure [2cl shows that, on average, >99% matching vox- 
els of a target heart object is reached after 6400 evaluations and >99.5% after 
24000 evaluations. 

At the end of the experiments, the fittest individual was subsequently fabri- 
cated by a three-dimensional printer and can be seen in Figure [3] including the 
supporting rafts required for manufacture. Figure [4] illustrates a sample of the 
evolved individuals from one cube experiment. Figure [S] similarly for the star 
experiment, and Figure [S] for the heart experiment. 

4 Conclusions 

This paper has shown that it is possible to evolve a vector of reals which are used 
as Superformula parameters to generate three-dimensional objects. Target-based 
evolution was used to explore the ability of Superformula to create complex ob- 
jects, particularly those that resemble natural designs. The experiments showed 
that with target-based evolution very closely matching objects can be identified. 
One significant advantage of the approach over alternative representations is the 
simplicity and compactness of the encoding. 
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(a) Cube mi = 4, ni_i = 
10, ni,2 = 10, ni,3 = 10, 
m2 = 4, ri2,i = 10, 712,2 = 
10, n2,3 = 10 




(c) Heart nii = 3, = 
1-5, ni,2 = 12, ni,3 = 3, 

m.2 = 0, 712,1 = 3, 712,2 = 0, 
712,3 = 




(e) Torus mi — 10, ni,i = 
10, ni,2 = 10, ni,3 = 10, 

7n2 = 10, 7l2,l = 10, 7l2,2 = 
10, 712,3 = 10 tl = 2, C3 = 




(b) Stax mi = 6, ni,i = 5, 

ni,2 = 10, 7ll,3 = 10, 7712 = 
4, 712,1 = 10, 712,2 = 10, 
7l2,3 = 10 




(d) Shell mi = 3, ni,i — 
1-5, ni,2 = 12, ni,3 = 3, 

mi = 0, 7121 = 3, 7122 = 0, 
7123 = 0, t2 = 2, dl = 1, 

d2 = 1, ci = 5 




(f) Mobius Strip mi = 3, 

711.1 = 1.5, 7ll,2 = 12, 
7ll,3 = 3, 7712 = 0, 7l2,l = 3, 

712.2 = 0, 7l2,3 = tl = 4, 

t2 = 0, rfi = 0, d2 = 0, 
ci = 5, C2 = 0.3, C3 = 2.2 



Fig. 1: Example three-dimensional Supershapes 
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Fig. 2: Evolution of three-dimensional Supershapes 




Fig. 3: Cube, Star and Heart fabricated by a three-dimensional printer. 



(a) 96.71% (b) 98.50% (c) 99.29% (d) 99.97% (e) Target 

Fig. 4: Evolution of a three-dimensional cube 




(a) 94.18% (b) 96.24% (c) 98.03% (d) 99.56% (e) Target 

Fig. 5: Evolution of a three-dimensional star 




(a) 96.11% (b) 98.11% (c) 99.00% (d) 99.50% (e) Target 

Fig. 6: Evolution of a three-dimensional heart 



